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Abstract. Mai'kov partitions work most efficiently for Anosov systems or for Axiom A systems. 
However, for hyperbolic dynamical systems which are either singular or whose hyperbolicity is nonuni- 
form, the construction of a Markov partition, which in these cases is necessarily countable, is a rather 
delicate issue even when such a construction exists. An additional problem is the use of a countable 
Markov partition for proving probabilistic statements. For a wide class of hyperbolic systems, L. S. 
Young, in 1998, constructed so called Markov towers, which she could apply successfully to estab- 
lish nice, for instance, exponential coiTelation decay, and, moreover, as a consequence, a central limit 
theorem. The aim of this survey is twofold. First we show how the Markov tower construction is ap- 
plicable for obtaining finer stochastic properties, like a local limit theorem of probability theory. Here 
the fundamental method is the study of the spectrum of the Fourier transform of the Perron-Frobenius 
operator These ideas and results are applicable to all systems Young has been considering. Second, 
we survey the problem of recurrence of the planar Lorentz process. As an application of the results 
from the first part, we obtain a dynamical proof of recurrence for the finite horizon case. Here basically 
different proofs were given by K. Schmidt, in 1998, and J. -P. Conze, in 1999. As another application 
we can also treat the infinite horizon case, where already the global limit theorem is absolutely novel. 
It is not a central one, the scaling is y'n log n in contrast to the classical ^/n one. Beyond thus giving a 
rigorous proof for earlier heuristic ideas of P. Bleher, which used three delicate and hard hypotheses, 
we can also a) verify the local version of this limit theorem for the free flight function and b) prove the 
recurrence of the planar Lorentz process in the infinite horizon case. 

1. Introduction 

Since — following some ideas of Hadamard — M. Morse introduced the concept of symbolic dy- 
namics, the method got more and more extensively used to study topological, and later also ergodic 
and stochastic, properties of dynamical systems possessing some hyperbolic behaviour 

On the one hand, "the idea of coding and semiconjugacies with topological Markov chains yields 
remarkably precise results concerning topological entropy, the growth of periodic orbits, the pres- 
ence of orbits of various periods, and the structure of maps with zero topological entropy" LKH 95 J . 
On the other hand, through the achievements of Bowen (cf. MB 751 . Ruelle (cf. BR 781 and of Sinai 
(most notably his work I'S 721 relating symbolic dynamics and Gibbs states of statistical physics), 
almost invertible semiconjugacies and conjugacies provided by Markov partitions made it possible 
to demonstrate exponential correlation decay and further nice and useful stochastic properties for 
Axiom A systems — and later for more general ones, too. Thus, for a long time it, quite naturally, 
seemed so that the construction of Markov partitions is 'the method' for obtaining effective statistical 
statements for more complicated systems as well. 
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However, in course of the work of the Moscow school on bilUards and of Benedicks and Young on 
the Henon map it became clear that a Markov partition is a too delicate construction if one wants to 
relax the assumptions of smoothness or of the uniform hyperbolicity of the maps in question. In fact, 
a warning might have come earlier from Bowen's result showing that, even in the nicest systems, in 
the multidimensional (d > 2) case the boundaries of the elements of any Markov partition behave 
wildly, in particular, they are not smooth fB 781. In addition, for instance, in discontinuous systems 
like billiards, the local invariant manifolds are, indeed, arbitrarily short, and as a consequence the 
Markov partition is necessarily countable and its elements are products of Cantor sets. Then to adapt 
the boundaries of a possible Markov partition in a Markov way is an extremely delicate issue even 
in the case when such a construction was successfully established (cf. I BS 801 ). 

It is not our aim to go into more details here since there exists a quite recent and excellent survey 
IIChY 01. which, on the one hand, gives a comprehensive historical overview, and, on the other hand, 
explains the way out: the 'weaker' construction of a Markov tower. This construction was designed 
by L. S. Young lY 981 and it works for a wide class of systems with some hyperbolicity, among 
others for Anosov and Axiom A systems, two dimensional hyperbolic systems with singularities, 
two dimensional Sinai billiards with a finite horizon, hyperbolic unimodal maps and hyperbolic 
Henon maps. 

The construction of a Markov partition and the resulting symbolic dynamics opened the way in 
a straightforward manner to put the probabilistic arsenal of — appropriately mixing — stationary 
stochastic processes into action. In the case of a Markov tower this connection is not straightfor- 
ward, and our actual aim is, indeed, to understand and to discuss the probabilistic approach in the 
case of a Markov tower. We note that, as we will see in subsection 2.2, the tower also leads to a 
(countable) Markov partition, but its properties, apart from its formal ones, are quite different from 
those of a traditional Markov partition. In particlular, it does not seem to provide a flexible Markov 
approximation. This is why in its applications new methods are needed and, in fact, their discussion 
is our main aim here. 

The paper is organized as follows: in section|2lwe recall the axioms of systems, which we are 
going to deal with, and briefly describe the tower construction. In section 3 we analyze how one 
can establish stochastic properties, in general, and further finer stochastic properties, like local limit 
theorems, in particular The choice of local limit theorems may seem eventual but it is not. This will 
be clear from section 4, where we apply the local limit theorems to planar dispersing billiards. In 
doing so, beside attaining local theorems for them, we also obtain 

• a dynamical proof of recurrence for the planar Lorentz process with a finite horizon (in 
fact, partly abstract ergodic-theoretic proofs were given by K. Schmidt liSch 98J and Conze 
ICon 991 . which were, however, also using the central limit theorem); 

• the first rigorous proof for a noncentral limit theorem for the displacements of the planar 
Lorentz process with an infinite horizon; this result was conjectured by an earlier nonrig- 
orous, heuristic argument of Bleher | B 92 1 based on three hard hypotheses (which, in fact, 
still do not follow from our approach); 

• the first proof of recurrence for the planar Lorentz process with an infinite horizon. 

2. Markov tower 

2.1. The Product Set. The technique developed in lY 981 allows to handle stochastic properties of 

systems 

(1) whose every power is ergodic; 

(2) which satisfy several technical assumptions well-known from hyperbolic theory; 
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(3) whose phase space X contains a subset A with a hyperbolic product structure; 

(4) where the return time into A has an exponentially decaying tail. 

This class contains planar dispersing billiards with both bounded or unbounded free flight (i. e. with 
a finite resp. an infinite horizon), logistic interval maps, expanding maps with neutral fixed points, 
piecewise hyperbolic maps, Henon attractors, their generalisations, and certain partially hyperbolic 
systems. The rest of this subsection is devoted to the precise definitions. 

We start with describing precisely the models we are going to deal with. Let T be a C^^*^ diffeo- 
morphism with singularities of a compact Riemannian manifold X with boundary. More precisely, 
there exists a finite or countably infinite number of pairwise disjoint open regions {X, } whose bound- 
aries are submanifolds of codimension 1, and finite volume such that UX, = X, T [^^^ is 1 — 1 and 
can be extended to a C'+^-diffeomorphism of Xi onto its image. The Riemannian measure will 
be denoted by jU, and if IV C X is a submanifold, then will denote the induced measure. The 
invariant Borel probability measure will be denoted by v. 

Definition 1. An embedded disk ycX is called an unstable manifold or an unstable disk ifix^y e 
Y, d{T^"x, T^"y) exponentially fast as n ^ oo; it is called a stable manifold or a stable disk 
ifyx,y G Y) d{T"x,T"y) — > exponentially fast as n^ oo. We say that F" = {y"} is a continuous 
family of C' unstable disks if the following hold: 

• K'^ is an arbitrary compact set; D" is the unit disk of some R"; 

• 0": xD" ^X isa map with the property that 

- <!>" maps x D" homeomorphically onto its image, 

- X <I>" I {{x} X D") is a continuous map from K" into the space ofC^ embeddings of 
D" into X, 

- y", the image of each {x} x D", is an unstable disk. 
Continuous families of C' stable disks are defined similarly. 

Definition 2. We say that AcX has a hyperbolic product structure if there exist a continuous family 
of unstable disks F" = {"f} and a continuous family of stable disks P — {y*} such that 

(i) dimY" + dim/ = dimX 

(ii) the -disks are transversal to the Y -disks with the angles between them bounded away from 
0; 

(iii) each y"-disk meets each Y-disk in exactly one point; 

(iv) A= (UY")n(Uf ). 

Definition 3. Suppose A has a hyperbolic product structure. Let F" and F* be the defining families 
for A. A subset Aq C A is called an s-subset i/ Aq also has a hyperbolic product structure and its 
defining families can be chosen to be F" and Fq with Fq C F''; M-subsets are defined analogously. 

For X € A, let"f{x) denote the element of T" containing x. 

Definition 4. We call (X, T, v) a Young system, if the following Properties (P1)-(P8) are true: 

(PI) There exists a A c X with a hyperboUc product structure and with /iylYfl A} > for every 
YeF". 

(P2) There is a countable number of disjoint s-subsets Ai , A2, • • • C A such that 

• oneachY"-disk^y/{(A\UA,)nY"} = 0; 

• for each i, 3Ri G Z+ such that r'^' A,- is a M-subset of A; 

• for each n there are at most finitely many i's with /?, = n; 

• min/?i > some Rq depending only on T 
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(P3) For every pairx,y e A, we have a notion of separation time denoted by so{x,y). If s(,{x^y) = 
n, then the orbits of x and y are thought of as being "indistinguishable" or "together" through 
their iterates, while r"+^x and T^^^y are thought of as having been "separated." (This 
could mean that the points have moved a certain distance apart, or have landed on opposite 
sides of a discontinuity manifold, or that their derivatives have ceased to be comparable.) 
We assume: 

(i) *o > and depends only on the y'-disks containing the two points; 

(ii) the number of "distinguishable" n-orbits starting from A is finite for each n\ 

(iii) forx,>'GAi, SQ{x,y)>Ri + so{T'^ix,T'^'y); 

(P4) Contraction along f disks. ThereexistC>Oanda< 1 such that foryGf(x), d{T"x,T"y) < 
Ca" Vn>0. 

(P5) Backward contraction and distorsion along 7". For y € 7"(x) and < A: < n < io(-^)>'). we 
have 

(a) d{T''x,T"y) < Ca*o(A:,)')-n. 
(b) 

(P6) Convergence of I7") and absolute continuity of P. 

(a) forj e f{x), 

IoeUT^VI<C^" Vn>0. 

(b) for y, y e r", if : yn A ^ y n A is defined by &{x) =Y{x)ni, then is absolutely 
continuous and 

^(Q.'W) / .fr detZ)P'(r'x) 
^A^Y i=o dstDT''{T'&x) ■ 

(P7) 3Co > and 60 < 1 such that for some y e T", 

ILiy{x e y n A : /?(jc) > n} < CoSg Vn > 0; 
(P8) {T",v) is ergodic Vn > 1, where V is the SRB-measure corresponding to the system {X, T). 

Remark 5. Some explanation: Properties (Pl-7) do not involve the measure V. Using them, in fact. 
Young constructs the SRB-measure v, so the system satisfying (Pl-7) is a Young sysstem if for the so 
constructed SRB-measure (P8) holds, too. 

2.2. The Tower. Now we wiU define the Markov extension, the actual Markov tower. Let /? : A — > 
Z+ be the function which is /?, on A,, and let 



and define 



A = {(x,/):xeA; / = 0,1,. . . - 1} 



^ \ {T\0) if l-\-l=R{x) 



We will refer to A; = {{x,j) \ {x,j) G A, j = 1} as the level of the tower A. Young also has a 
construction for v, the SRB measure of the extension, for which the pushforward is V, and J{F) = 1 
except on (Aq). Thus, the Markov Tower is the dynamical system (A,F, v). 
On the tower a Markov partition D can be defined, with the following properties: 
(a) D is a refinement of the partition A;. 
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(b) By denoting by ©/ the partition ® |A/, 2)/ ~ {Aij \ j — I, . . . has only a finite number of 
elements and each one is the union of a collection of A,'s; 

(c) 2)/ is a refinement of 

(d) if X and y belong to the same element of 2)/, then SQ{F^'x,F^'y) > I; 

(e) if Ri = Rj for some / j, then A, and A^ belong to different elements of 'Dr.^i. 

Let A'l J = A/.j nF"'(Ao). We think of Aij\AJ j as "moving upward" under F, while A^ j returns to 
the base. 

def 

It is natural to redefine the separation time to be s{x,y) ~ the largest n such that for all / < «, F'x 
and F'y lie in the same element of {A/ j}. We claim that (P5) is valid for x,y £ 7" fl A/ j with s in the 
place of SQ. To verify this, first considers, y € A. We claim that s{x,y) < so(x,y). If x,y do not belong 
to the same A,-, then this follows from rule (d) in the construction of 2)/; if x,y G A,-, but T^x, T^y 
are not contained in the same Ay, then s{x,y) = Ri + s{T'^x,T^y), which is < so{x,y) by property 
(P3),(iii) of 50, and so on. In general, forx,y e A;./, let xq — F^'x, yo = F^'y be the unique inverse 
images of x and y in Aq. Then by definition s{x,y) ~ s{xo,yo) — I, and what is said earlier on about 
XQ and yo is equally valid for x and y. 

From here on sq is replaced by s and (P5) is modified accordingly. 

3. Stochastical Properties 

As explained above, the tower construction actually provides a countable Markov partition. Be- 
low we first remind the reader how Young exploits this partition to obtain stochastic behaviour The 
first goal is to estabhsh exponential decay of correlations, where, of course. Property (P7) plays a 
decisive role. 

3.1. The Perron-Frobenius Operator and the DoebUn-Fortet Property. 

1. The starting point is to facto rise the dynamics by a factorisation along stable manifolds of A. The 
advantage is that this dynamics will behave as an expanding map, an appropriate object to study via 
the Perron-Frobenius operator. Let A := A/ ^ where jc ~ y iff y S Yi^)- Since F takes "/-leaves to 
y^-leaves, the quotient dynamical system F : A ^ A is clearly well defined. 

The construction of fh is not trivial but is quite standard. Young obtains it following IB 751 . A 
simple property of m is: let ;n|A/ be the measure induced from the natural identification of A/ with a 
subset of Ao, so that /(F) = 1 except on F'^Aq), where J{F) = 7(7^ o We note that for 

the factorised map Young also proves a distorsion property with a weaker constant P : 1 > p > 

To investigate Birkhoff sums we have to associate a function / : A ^ M to each observable / : 
M ^ M. We can pull back / to the tower, and find an other function cohomologous to this one, 
which is constant along stable manifolds. This method is described for example in 1 PP 90J . 

2. The analytic tool of investigation is the Perron-Frobenius operator: 

y:Fy=x ' 

3. The technique is based upon the spectral properties of the Perron-Frobenius operator. For this 
purpose we have to introduce suitable Banach spaces, where P has a nice spectrum. Actually for the 
method introduced by Doeblin and Fortet we need a pair of function spaces C and L (usually with 
some supremum-like and Lipschitz-like norms) such that iS<C,|l-|lc<|llU, and the inclusion of 
L into C be a compact operator 

If we have such a pair of Banach spaces, and we can prove, that 3N,K, and T < 1 



l^^'^IU <^ll9lU+^ll9ll 



c 
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then by knowing that V/ T' is ergodic we have that the spectrum of P on L is contained in a disk 
with radius strictly smaller than one, except that 1 is an eigenvalue with multiplicity one, and the 
corresponding eigenfunction is the invariant density (cf. [I-TM 501). 

This kind of estimate captures the uniform expanding feature of the dynamics, or contraction of 
the P operator. In order to derive this so called Doeblin-Fortet property of the transfer operator (in 
the theory of dynamical systems often called the Lasota-Yorke property, cf. fLY 731), Young uses 
an exponential factor in the function norms. F on the tower has Jacobian 1, when moving upwards, 
and the tower is usually unboundedly high, so we have to pretend expanding at least in the norms: 



The denominator in the second definition is a natural distance for the points on the tower, so it 
is really a Lipschitz-like norm. The aforementioned trick is in the e^'^ term. This allows that a 
function on A which is exponentially increasing with the height of the tower, to be in C, if this 
growth is moderate. This also means, that P is contracting the norm in the middle of the tower, 
where the Jacobian is 1 . The constant e should be chosen carefully in order to hold back enough 
contraction when a Markov return occurs. Roughly speaking £ has to be smaller than the smallest 
positive Lyapunov exponent. 

3.2. The Central Limit Theorem. The aforementioned spectral picture is essentially equivalent to 
the exponential decay of correlations for function pairs in X . As a matter of fact, the boundedness 
of the functions in question is also needed. After the argument presented in | Y 98 1 this leads to the 
exponential correlation decay of bounded, piecewise Holder functions on M. For the same class of 
functions Young immediately gets the Central Limit Theorem (CLT) by checking the conditions of 
a theorem by Keller |K 80 1. Here we make a simple but important clarification. The traditional and 
by far the most widely used method for establishing the CLT is to apply Fourier transforms. Keller 
||K 80 1 (and thus also | Y 98|) can elude this by referring to a nice and useful theorem of Gordin 
llG^69i . (As a matter of fact, Fourier transform are, indirectly, still applied, since Gordin constructed 
a martingale approximation and used the martingale CLT. But in proving the CLT for martingales 
again Fourier transform is the method! Moreover, the error term in the martingale approximation is 
so large that it obviously excludes the applicability of this approach in proofs of finer statements, for 
instance, in those of local limit theorems.) 

Let we recall two basic results from I Y 981 : the first one on the exponential decay of correlations 
and the second one on the CLT. Notation: 

i?fTi = {(p :M^M I 3A > OsuchthatforVx,y eM |(p(jc)-(p(y)| <Ad{x,y)'^] . 

Theorem 6. f lY 981 ) For any r|, there exists T < 1 such that for all (p,\|/ G there exists a C = 
C((p,\(/) such that 



sup (p 



e 





.n 



Vn > 1 



Theorem 7. ( lY 981 ) Assume (p G i?/n and J cpc/v = 0. Then 
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for some O > 0. Moreover, O — Oiff<^ = \\foT — ^ for some \(/ G i"(v). (Here =^ denotes weak 
convergence of probability distributions and 0\i {m,0^) denotes the normal distribution with mean m 
and variance O^.) 

For simplicity, we have formulated these results for one-dimensional random variables, and their 
extension to vector valued functions is straightforward. 

3.3. Local limit theorem. For illustrating a local CLT as compared to the widely used (global) CLT, 
take a simple symmetric random walk (SSRW) on l/. So let W„ = Xi + . . .X„, where Xi,. . . ,X„, . . . 
are independent, identically distributed random variables with the common distribution P{Xi = 
±ey) = i < j < d for all / e Z+ (here the e,s are the standard unit vectors of Z''). Then, of 
course, the CLT says that P{W„ e ^/nA) ^{A) as n ^ °o, where 0(A) — J^<^{s)ds and (^{s) = 
(2jt)^''/^exp(— y) is the li-dimensional Gaussian density. In other words it describes the asymp- 
totics of a sequence of sets increasing like y/n. In contrast, for the SSRW the local CLT (LCLT) says 
that n''^^P{W„ — [s^/n\) — > (^{s) as n ^ i. e. it describes the asymptotics of a sequence of sets of 
fixed size, consequently it is, indeed, local! 

For stating our main theorem we have to fix some notations first. For a fixed f : X denote 
the average / fdv — a, and 

S{ix)^"ff{T'^x) 

k=0 

the Birkhoff sum. Consider the smallest translated closed subgroup V + r CM.'' which supports the 
values of / {V is the group and r is the translation). By ergodicity of all powers of T, the support of 
S„ is y + nr. 

Theorem 8. f' ISzV 031 Suppose that 

(1) {X, r,v) is a Young system (cf. subsection \2.1\ : 

(2) / is minimal: i. e. it is not cohomologous to a function for which the support in the above 
sense is strictly smaller 

(3) / is nondegenerate: i. e. span (V) — W', and 

(4) / is bounded and Holder continuous. 

Let k„ + nr be such that k. Denote the distribution of St — k„ by X)„. Then 

n''l^-V„ -^(^{k)l 

where (|) is a nondegenerate normal density function with zero expectation, and I is the uniform 
measure on V: product of counting measures and Lebesgue measures. The convergence is meant in 
the weak topology. 

Remark 9. For nonminimal functions we can obtain an analogous result. The limit measure on the 
right hand side in this case is not necessarily uniform. 

We know from classical analysis or from probability theory that the local behaviour of distri- 
butions (densities, measures,...) is connected to the tail behaviour of the corresponding Fourier 
transforms. Therefore, in the next subsection, we are going to study the Fourier transform of the 
Perron-Frobenius operator. 
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3.4. The Fourier Transform of the Perron-Frobenius Operator. When one wants to obtain finer 
results tiian the CLT, then Fourier transforms seem inevitable. Thus, in our setup, we have to define 
the Fourier transform of the Perron-Frobenius operator: 

where / is the function for which the limit theorem is stated. Note that Pq = PI Also, it is worth 
noting that P/'l = P"(exp(/f5,()) and, moreover, EP,"1 = EP"(exp(/r5,()) = exp{itsl,), the usual 
characteristic function of sl where s{, denotes the Birkhoff sum for the function /. 

The heart of this method is to expand the leading eigenvalue of P, — analogously to the Taylor 
expansions around of characterisctic functions of probability theory. Before that, however, we 
have to ensure its existence. We proved in |SzV 03 1 that 1 1-^ Pt is continuous in the il-norm, and 
by the stability of the spectrum, for small values of t there exists A-r, the perturbed value of 1 as an 
eigenvalue with multiplicity one. By proving the Doeblin-Fortet inequality for P, we get that the 
rest of the spectrum will lie in a disk, with radius smaller than 1, so it will not bother to be the 
leading eigenvalue. 

If the function / is bounded and piecewise Holder continuous on M, then we get the second order 
Taylor expansion for the Fourier transform: 

P,((p) =P((P) +/fP(/(p) - ^P(/-2(p) +o(f2) ||/2(p||^ . 

From the assumptions it follows that is finite. By an argument presented in various forms 

in |,Nag 57 J , LKSz 831 and [GH 88 1 this leads to the expansion of A,,. 

The philosophy explained above can already be combined with the classical proof of the local 
CLT sketched in the Appendix. Indeed, term // is the same, term / corresponds to the CLT. Term /// 
can be handled by using the ideas outlined above, while for bounding term //// we have used some 
compactness arguments borrowed from IAD Oil . 

4. The Planar Lorentz Process 

4. 1 . Semidispersing billiards and Lorentz process. In this subsection we summarize some basic 
properties of semidispersing billiards. Our aim is to introduce the most important concepts and fix 
the notation which is essentially borrowed from IKSSz 901 . Semidispersing billiards are more or 
less hyperbolic dynamical systems with singularities. Pesin's theory was extended to these systems 
in IKS 861 . 

A billiard is a dynamical system describing the motion of a point particle in a connected, compact 
domain Q cT'' . The boundary of the domain in assumed to be piecewise -smooth. Inside Q the 
motion is uniform while the reflection at the boundary dQ is elastic. As the absolute value of the 
velocity is a first integral of motion, the phase space of the billiard flow is fixed as M = Q x S''^ ' 
- in other words, every phase point x is of the form x — {q,v) with q E Q and v G M'', |v| = 1. 
The Liouville probability measure /j on M is essentially the product of the Lebesgue measures, i. e. 
d/i = const. dqdv. The resulting dynamical system {M,S^,ij) is the (toric) billiard flow. 

Let n{q) denote the unit normal vector of a smooth component of the boundary dQ at the point 
q, directed inwards Q. Throughout the paper we restrict our attention to semidispersing billiards: 
we require for every q e dQ the second fundamental form K{q) of the boundary component to be 
nonnegative. 

The boundary dQ defines a natural cross section for the billiard flow. Namely consider 

dM ^ {{q,v) \q edQ,{vMq)) >0}. 
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This set actually has a natural bundle structure (cf. IBChSzTI '). The Poincare section map T, also 
called the billiard map is defined as the first return map on dM. The invariant measure for the map 
is denoted by /Ji, and we have dp\ = const. | {v,n{q))\dqdv. Throughout the paper we work with this 
discrete time dynamical system {dM,T,/Ji). Recall the usual notation: for {q,v) e M one denotes 
n{q, v) ^q the natural projection. 

The Lorentz process is the natural 7/' cover of a toric billiard. More precisely: consider n : 
R"^ — > T"^ the factorisation by Z''. Its fundamental domain D is a t/-dimensional cube (semi open, 
semiclosed) in R'', so R'' — Ll,^^d (D + z), where D + zis the translated fundamental domain. 

By denoting Q = U'^Q, M = Qx S''^^, etc., the Lorentz dynamics is {M,{S' \ t € R},ft) and 
its Poincare section map is (dM,f The free flight function \\i : dM W' is defined as follows: 
Tij/(i) = q{Tx) — q{x). The discrete free flight function ic : dM — > Z'' is defined as follows: k{x) = 
\{tx) — where i(x) = z if i G Dz- Observe finally, that \j/ and ic are invariant under the l/ 
action, so there are \(/ and K functions defined on dM, such that x}/^ = n*\(/ and K = n*K. Actually 
for our purposes it will be more convenient to choose the fundamental domain in such a way that 
dQ{^dD = 0. In this way K will be continuous. 

4.2. LCLT and recurrence for the d — 2 finite liorizon case. Consider the Lorentz process start- 
ing from the fundamental cell, a fixed isomorphic version of the fundamental domain. In the domain 
the starting phase point is random, it is distributed according to the invariant measure of the corre- 
sponding torus-billiard. The relative position of the Lorentz particle after the «* collision is St, the 
discrete position is 5]^. So the event A„ = (5,^ = 0) means that after the n* collision the particle is 
again in the fundamental domain. Recurrence means that it happens almost surely. 

Tlieoreni 10. ISch 981 . BCon 991 . BSzV 031 The planar Lorentz process with a flnite horizon is 
recurrent. 

For the proof we will use a stronger version of the well-known Borel-Cantelli lemma. This 
version is due to Lamperti: 

Lemma 11. f USpi 64| j If for the sequence of events {A,,} 

£ v(Ai) = oo 

k=\ 

and some asymptotic independence holds: 

.1 v(A,A,) 

liminf — j < c 

(e v(...,) 

then An happens inflnite often almost surely: 

V^(^,v)e8M Bn^^oo (^,y) e p|A„,j = I 

To check the first condition we have to deal with the asymptotic probabilities of the events Sf^ = 0. 
This is exactly the region of the local limit theorem. Restrict ourselves to the d — 2 finite horizon 
case! This case is known to be a Young system i. e. it satisfies (P1)-(P8), so the first condition of 
our LCLT is satisfied. For the second condition we need to check the minimality of k: 
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Theorem 12. f ISzV 031 ) k is minimal in the class o/v|/. 

The proof is quite involved. Surprisingly it is related to arguments in IBChS 91 1. fBSp 961 and 
|B_00 I for establishing the nondegeneracy of the covariance matrix in the CLT To prove minimality, 
for each sublattice L C of finite index n = L : I?, we needed a periodic point n|per(.Ji:) such that 
the Birkhoff surn5^^^;c) ^ L. To find this point we used again Markov properties of A. Details 
can be found in ISzV 03l . 

Since K is minimal and the values span the plane it is nondegenerate. Since it is piecewise constant 
it is also Holder continuous, so the LCLT applies. This gives v{S^ = 0) ^ and the sum clearly 
diverges. 

The second condition of Lamperti's Borel-Cantelli lemma contains intersection of recurrence 
events, so we also had to prove a LCLT for joint distributions. 



5. Infinite Horizon 

Chernov ICh 991 observed that the planar infinite horizon Sinai billiard also satisfies (P1)-(P8), 
so according to theorem |8] the LCLT also holds exactly the same manner as stated for the finite 
horizon case. Nota bene: for bounded functions f : M —^W' satisfying the conditions in Theorem|8] 
the asymptotics of the Birkhoff sum v{sl = 0) ~ -jjj- 

But the \|/ or the K free flight functions are not covered by this theorem since here K is not bounded, 
\|/ is even not Holder This is not a surprise, since in this case a new phenomenon shows up. For- 
mer heuristic arguments by Bleher |B 92| already suggested to expect that the moving particle is 
superdiffusive namely ^^^'^^^ will have a limit distribution. The reason is the following: 

The infinite horizon condition is equivalent to the existence of collision free orbits in the phase 
space M. These orbits form corridors in R^. Large free flights can occur by "crossing" one of these 
corridors. The smaller the angle with the direction of the corridor is, the longer is the free flight. 



Computing the invariant measure for the small angle sets gives the asymptotics v(|k| = u) ~ ^. 
This means that / |Kpiiv = °°, but any power with smaller exponent is integrable. As a matter of 
fact, the distribution of K is in the nonnormal domain of attraction of the normal law, and its Fourier 
transform is: 

K = j e''Vv = l+c|f|2log|f|+0(|f|2) 

(where c is a constant matrix) which means that if we would add independent copies of the same 
distribution then , ," would tend to a gaussian law. The log factor comes from the fact that the 

\/nlogn ° ° 

truncated k'' = k1|k:|<x has a variance of order logx. 

Direct geometrical calculations show that, when K is large, then the order of K o r is between \/k 
and K^. (There is a case when the trajectory hits the neighbouring scatterer on the same side of the 
corridor before "crossing". 
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In this case we change the Poincare section and consider the sum of the small and the necessarily 
large free flight vector as KoT.) An important observation is that typically the next free flight will 
be in the regime of \/k- More precisely, for any 5 > 5 the probability v(ko T > | k = m) as 
u °o. Even more the conditional expectation Ev(ko r | k) is of order This means that though 
|Kp is not integrable, nevertheless, / k(ko T)dv < °°. The hope that the autocorrelation may have 
a fast decay has lead to the conjecture that asymptotically behaves the same way as the sum of 
independent copies of K. 

5.1. Limit theorems: global and local, and recurrence for the d — 2 infinite horizon case. To 

reach the aforementioned limit theorem, and moreover the local limit theorem for K we used the 
symbolic space A as constructed in IIY 981 . and investigated carefully how and where large values 
of K appear on the tower Clearly K is bounded on A. Since in one step it can grow at most to it is 
square, or shrink to it is square root the time needed to reach the K = m set from A is about c log log u, 
and before returning to A also the same amount of time is needed. By (P7) it is immediate that both 
the measure of phase points, which spend k iterates in the corridor before returning to the base of 
the tower, and both the measure of phase points visiting the corridor k times before returning to the 
base is exponentially small in k. 

We had to replace the function norms in the definition of C and L to refer to ic on the tower 
Instead of Young's e^^' factor we used 

f[min(e-^K-^oF-M . 

k=i ^ ' 

Observe, that if K is bounded and 5 is small enough this gives back e^^'. In that way we managed 
to achieve that f P; be a countinuous mapping in both function norms. This was not the case with 
the original norms, if the horizon is infinite. So P, can be considered as a perturbation of P, and this 
also extends to the leading eigenvalue A,,. 

To obtain that the asymptotic behaviour is the same as in the case of the sum of independent 
copies, we needed to show that 

X, = l+c|f|2log|f| + 0(|f|2) 

has the same kind of behaviour as the dynamically untouched ic. This is the key of the proof. While 
sums of independent copies give the product Fourier transform ic", that of the Birkhoff sums give 
/ P"(l)iiv. The operator can be approximated by A-"v (here v is the leading projection operator) up 
to an exponentially small error term coming from the rest of the spectrum. Summarising: if kit) 
and X, behave the same way as f ^ 0, then the sum of independent copies and the Birkhoff sum Sf^ 
behave the same way as n °o. 

There is a quite involved proof of the A,, expansion which is based on correlation estimates of 
powers of the truncated k and the eigenfunction related to A-;. Details will appear in a technical 
paper. 
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Theorem 13. Suppose that the direction vectors of infinite collision free flights span the plane. Let 
A C M^. Then 

y{S^ e JnlognA) ( (Sf{k)dk 

J A 

where (|) is a nondegenerate gaussian density with zero expectation. 

Remark 14. The problem of the limiting behaviour of displacements in the case of an infinite horizon 
has raised the interest of several people using different methods (very interesting works are IB 921 
and |ZE 91 \). It is worth mentioning that the computational method of l>ZE 97 J forecasts a non- 
Gaussian limit under the same scaling. . 

Theorem 15. Suppose that the direction vectors of infinite collision free fiights span the plane. Let 
k„ G such that , \" — > k. Then 

nlog« v{Sl = k„)^m 
where (|) is a nondegenerate gaussian density with zero expectation. 

isf depends only on the corridor geometry. When computing the covariance one considers only 
the directions and widths of corridors, and the bounding points of the corridors. This is a finite set 
of points on the scatterers, the geometry of this finite set, and the curvature of the scatterers at these 
points are involved, but nothing else. The sum of -^^^^^ diverges so the recurrence follows using an 
analogous argument as for the finite horizon case. In the case when all corridors are parallel one has 
to apply a nonisotropic scaling. 

Theorem 16. Suppose, that all collision-free flights in the plane are parallel to the unit vector w. 
Consider the linear transformation Bn which has the matrix ^ ^"q°^" ^ ^ in the basis w,w-^. Let 
k„ S 1? such that B^^k„ k. Then 

detB„ viS"^ = k„) = n^/l^ v{Sf, = k„) ^ ^ 
where (|) is a nondegenerate gaussian density with zero expectation. 

The sum of — J — is also divergent, so recurrence is also obtained in this case. 

Remark 17. It is worth noting that in a recent manuscript of Gouezel IG 021 a related problem 
was investigated for \ — D piecewise expanding maps with a neutral flxed point. The behaviour of 
billiard orbits in corridors is analogous to that of orbits of these \ —D maps near the neutral flxed 
point. Essential difficulties in our case arise from a) the larger dimension of the space; b) the not 
quite explicit form of the billiard dynamics; and c) emphatically from the fact that in our case the 
function of interest is unbounded with a quite long tail whereas in IG 021 it is bounded. His setup is 
more general since he is also considering stable limit laws in general, like IAD Oil . The restriction 
of our interest to the nonnormal domain of attraction of the Gaussian law came from the fact that 
our main concern was the free flight function. 

Acknowledgement. The authors thank the refereee for his valuable remarks and suggestions. 
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Appendix: A Classical Local CLT 



Here we recall, in the simplest case, the classical proof of Gnedenko IG 481 of a local CLT. 
For simplicity, consider the case d = I. Following the notations of subsection l3.3l our goal here 
is to prove that, as n oo, 

if ^ k. Heuristically one expects that 

VnPiWn = kn) = VnP{^ = -^) ~ = - / cy^pi- is ^)j{s)ds 

2 

where y(s) = exp(— y) is the standard gaussian characteristic function. Here we used the gaussian 
approximation suggested by the CLT and the Fourier inversion formula. 

For a proof, let us turn to characteristic functions of W„ . Denote by ^ (f ) the common characteristic 
function of the variables X,. Then 



and by substituting s = t^/n, this is equal to 



- 1 

I — 

271 



exp(-/f^„)^"(f)(if 



1 

27t 



exp(/s-^)^"(-^)c/i 



We emphasize that it is fundamental to precisely know the support of the values of the variables X, 
(i. e. the minimal lattice containing these values) since the form of the inversion formula depends 
on this. Moreover, the fact that the X,s take their values on a lattice of span 2 was used in our first 
heuristic formula, too. (This information is encapsulated in the minimality condition of Theorem 7.) 
To prove our desired statement we write 



< 



s\<A 



ds - 



2% 



s>A 



V^P{Wn = k, 

'Y{s)ds 
= 1 



-2^(- 



A<|.si<ev^ 

II + III + nil 



ds- 



ds 



For making the right hand side sufficiently small, we will first select A to be sufficiently large and 
then e sufficiently small. Thus // can be made arbitrarily small, and, for fixed A, / will also be small 
by the CLT (as a matter of fact, the smallness of / would also follow from our forthcoming argument 
for handling ///. In fact, /// is a quite interesting term. By expanding t,{t) in a power series in the 
neighbourhood of 0, one can easily see that, if e is sufficiently small, then for |s| < e one has 



I^"(^)I<1- 



4 <exp(--) 



As a consequence one obtains 



UK 



_exp(- 



-)ds< [ 
Ja< 



exp( )ds 
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and this term is also small if A is large. Finally, the smallness of //// follows from the fact that, in 
the interval e^/n < \s\ < the term ^(■^) is uniformly bounded away from 1 from above, and 
consequently for the integrand in //// we have an exponentially collapsing upper bound. 
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